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ABSTIU\CT 

We  consider  a  gas  composed  of  atoms  having  two  electronic 
states.  Transitions  between  these  two  states  ccn  take  place  by  collis- 
ions, and  by  emission  or  absorption  of  radiation.  When  the  radiative 
transition  probabilities  are  small  compared  with  those  due  to  collisions, 
a  •  temperature"  may  be  defined  for  the  translational  motion.  We  formu- 
late the  general  problem  for  the  case  where  the  total  concentration  of 
the  gas  depends  on  the  space  co-ordinates  in  a  (p.ven  w&y,  and  the  possibi- 
lity of  the  "  imprisonment"  of  the  radiation  emitted  is  taken  into  account, 
In  "this  case  the  problem  leads  to  a  pair  of  coupled  integro-differential 
equations,  from  which  the  "  temperature"  and  the  concentration  of  atoms 
in  the  excited  state  are  obtained  as  functions  of  space  and  time.  Neglect- 
ing the  resonance  radiation  and  the  spatial  variation  of  the  temperature, 
we  have  made  calculations  for  the  case  when  the  energy  difference  between 
the  two  states  is  of  the  order  kT,  for  a  number  of  values  of  the  ratio  be- 
tween radiative  and  collisional  transition  probabilities.  The  results  are 
applied  to  the  problem  of  the  cooling  of  the  atmospheric  gas  at  high  altit- 
udes (100  Km)  at  night,  as  a  result  of  the  magnetic  dipole  transitions  be- 

3 
tween  the  components  of  the  P  state  of  the  oxygen  atom. 
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1,  Introduction 

In  a  recent  paper,  Bates  -'  sugf^ested  thst  magnetic  dipole  transi- 

3 
tions  between  the  components  of  the   P  state  of  the  oxygen  atom  could  lead 

to  a  significant  coolinf^  of  the  upper  atmosphere  in  that  region  where  oxygen 
atoms  are  abundant,  i.e.,  about  100  Km  above  the  earth.  The  (inverted)  tri- 
plet state  in  question  has  splittings  of  the  order  of  ,02  e,v,,  which  are  of 
the  order  of  the  gas  kinetic  energy  at  these  altitudes,  so  that  the  number 

of  atoms  in  excited  states  of  the  triplet  will  be  largej  and  thus,  while  the 

3     3     3     3 
transition  probabilities  of  the  P_^  -   P,  ,   P,  -  Pp  transitions  are 

small  (  10~  -  10~  sec"  }   see  Pasternack'-  -'),  yet  the  total  number  of  transi- 
tions per  unit  volume  per  second  may  be  large  and  so  may  be  the  total  energy 
transferred  in  this  way. 

These  considerations  are  limited  by  the  ;-ssumption  of  a  Boltzmann 
distribution,  which  is  valid  only  for  a  system  in  thermodynamic  equilibrium, 
and  by  the  neglect  of  the  possibility  that  the  radiation  emitted  will  be  ab- 
sorbed by  other  atoms,  leading  to  the  familiar  problem  of  the  ••  imprisonment 
of  resonance  radiation"  ,  The  purpose  of  the  present  work  is  first,  to  formu- 
late the  problem  of  the  cooling  of  a  gas  by  the  emission  of  radiation  without 
assuming  the  Boltzmann  distribution  for  the  atoms  in  various  states  and  without 
neglecting  the  effect  of  imprisonment  of  the  radiation,  and  second  to  carry  out 
numerical  calculations  appropriate  for  application  to  the  upper  atmosphere  when 
the  effect  of  imprisonment  of  the  radiation  is  neglected  as  a  first  approxima- 
tion. 

To  simplify  the  calculations,  we  shall  consider  a  gas  composed  of  atoms 
having  two  electronic  states  1  and  2  ,  and  distributed  in  space  in  a  given 
manner  (such  as  the  gas  in  the  atmosphere).  Transitions  between  these  two  states 
can  take  place  by  collisions  and  by  emission  or  absorption  of  radiation.  In  a 
collisional  transition,  energy  is  transferred  between  the  translational  snd 
electronic  modes  of  the  atomj  in  a  radiative  transition  an  atom  may  emit  a 
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quantum  hv ^  of  radiation,  or  it  may  absorb  a  qusntun  tft/^^  that  has  been 
emitted  by  another  atom  in  the  gas.  We  shall  assume  that  there  is  no  incident 
radiation  from  an  external  source  having  a  frequency  v^  ,  so  that  the  radia- 
tion energy  of  this  frequency  comes  only  from  quanta  emitted  by  atoms  in  the 
system.  These  quanta  will  eventually  be  lost  to  the  system  after  undergoing 
a  series  of  absorptions  and  re-emissions  by  other  atoms  in  traversing  the  gas. 
The  problem  now  is  to  study  the  change  of  state  of  the  system,  starting  with 
suitable  initial  conditions.  At  time  t  =  0  our  gas  is  in  statistical  equili- 
brium at  a  given  temperature  T   as  a  result  of  collisions  alone.  Then,  for 

t  >  0,  we  "  switch  on"  the  radiative  transitions,  so  that  there  is  a  net  loss 
of  energy  from  the  system  for  t  >  0  ,  and  we  now  ask  the  following  question: 
at  whet  rate  is  energy  lost  from  the  system? 

The  behavior  of  the  system  obviously  depends  on  the  spatial  distribu- 
tion of  the  atoms,  and  on  the  values  of  the  collisional  end  radiative  cross- 
sections.  Let  the  average  cross-section  for  gas  kinetic  collisions.be  (JT,  and 
let  those  for  excitation  and  de-excitation  of  the  state  2  by  collisions  be 
C"    and  (Tp^  ,  respectively,  and  let  the  Einstein   A  coefficient  for  the 
2-1  transition  be  A^^  ,  Let  the  total  number  density  of  the  gas  be  N 
and  the  average  velocity  of  the  atoms  be  v  ,  In  the  extreme  case  where 

*21  ^  ^*^'  ^21  ^-"-^ 

the  excited  atoms  in  state  2  are  removed  by  emission  of  radiation  as  soon 
as  they  are  excited  by  collisions.  The  rate  of  loss  of  energy  of  the  system 
is  governed  by  the  rate  of  collisional  excitation  (which  is  a  function  of 
time),  and  by  the  effect  of  the  "  imprisonment"  of  the  radiation  in  the  gas. 
In  this  case,  the  distribution  of  the  atoms  in  the  two  states   1   and   2  , 
which  varies  with  time,  will  not  be  given  by  the  Boltzmann  theorem  at  all  since 
the  departure  from  thermodynamic  equilibrium  is  so  great  that  it  is  not  meaningful 
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to  speak  of  a  "  temperature"  .  In  the  other  extreme,  where 

Ag]^  «  N.v.  Cr^]^   ,  (2) 

the  rate  of  loss  of  energy  by  radiation  is  so   small  that  one  may  regard  the 
system  as  undergoing  a  continuous  transition  through  states  which  differ  only 
slightly  from  the  equilibrium  state. 

These  extreme  cases  can  be  treated  easily,  but  are  not  general  enough 
for  our  purpose.     In  the  present  vrork  we  shall  consider  the  ccse  where 

^21  ^^  ^'^^  ^o  '  ^^^ 

which  is  less     stringent  than  (2),   since  in  general    G"    >  "Jp-,    >  and  usually 
C  »  (Tot    •     Under  the  condition  (3)>   the  frecuency  of  elastic  col].isions 

O     ii-i- 

is  high  compared  with  that  of  emission  of  radiation,  so  that  the  atoms  make  enough 
collisions  to  redistribute  their  kinetic  energy  of  translation  between  radiative 
transitions.  In  this  case,  it  is  possible  to  assume  that  while  the  distribution 
of  the  atoms  in  the  two  electronic  states  1  and  2  may  deviate  greatly  from 
the  Boltzmann  distribution,  yet  the  translational  degrees  of  freedom  may  be 
characterized  by  a  "  temperature"  T  ,  defined  by  the  kinetic  energy  of  the 
atoms.  On  this  assumption,  the  problem  will  be  to  find  the  "  temperature"  T 
and  the  ratio  N^/N^   of  the  numbers  of  atoms  in  the  two  states   2  and   1 
as  functions  of  time  and  space  for  certain  initial  and  boundary  conditions. 

2.  Fonrtulation  of  the  genersl  problem 

We  consider  a  !7as  species  of  number  density  N  »  N(r)  that  has  two 
states,  a  ground  strte  1  and  an  excited  state  2  ,  of  densities  N.(r  ,  t)  , 
(i  =  1,  2): 


N(r)  -  N^(r,  t)  +  N2(r  ,  t)   .  (U) 
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We  shall  define  an  effective  temperature  T  ■  T(r,  t)  in  terms  of  the  total 
kinetic  energy  per  unit  volume.  It  is  convenient  to  deal  with  the  case 
i^en  there  are  other  gas  atoms  present  in  addition  to  the  gas  species  we  are 
concerned  with  primarilyj  these  other  atoms  enter  into  the  thermal  collisions, 
but  do  not  participate  in  any  way  in  the  radiative  processes.  (In  the  problem 
under  consideration,  there  would  be  nitrogen  molecules  present  as  well  as 
oxygen  atoms).  Accordingly,  the  total  number  density  of  all  particles  will 
be  X(r)  ,  which  may  be  greater  than  N(r).  Thus  the  total  kinetic  energy  per 
unit  volume  is 

(3/2)kTX  (5) 

which  serves  to  define  the  effective  temperature  T  • 

Let  the  col] isional  cross-sections  for  the  excitation  and  de- 
excitation  of  the  state  2  be  ^--(v)  and  (Tp-jCv),  where  v  is  the  relative 
velocity  of  the  colliding  particles.  There  will  also  be  radiative  effects 

specified  by  the  rates   B^g  ^^^2'  -'^^   ^^     {^21  *   ^2lP^'^21'  -'^A     ' 
where   A  ,  B   are  Einstein  coefficients  and  p  is  the  radiation  energy  density, 
2.1  Co!lJisional  transitions 

In  a  state  of  equilibrium  determined  by  collisional  processes  alone, 
we  have 

N.  =  N°(r)  only,  T  =  T°(r)  only,         (6) 

and 

N°(r)A°(r)  =  (g^/g-^)   exp  I  -   e^^Ar^C r)  |    (u) 

where  the  level  splitting   e^  is  given  by 

^21-  1^^^  -^^21  •  ^^^ 

IL.   is  the  reduced  mass  of  the  colliding  particles,  and  g,,  g„  are  the  statistical 
weights  of  the  states  1  and  2  . 
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If  X  is  the  total  number  density  of  particles  that  may  excite 
or  de-excite  our  atom,  the  rate  of  increase  of  N-  due  to  collisions  is 

CO  oo 

(  §1   N^)  ^^^  -  N^X  y   dv.v.f(v).  a--^^{v)   -  N^X  y  dv.v.f(v).  (J ^^{v) ,        (9) 

V  o 

o 

where  f(v)  is  the  velocity  distribution  function. 

The   two  cross-sections    C-in   t       Co-i    are  not  independent,  but  are 

related  atomic  properties.     In  particular,   if  we  define  two  velocities 

V  ,  v'    ,   related  by 

v^  =  v'^  +  (2/M^)   ep3_  =  v'^  +  V  ^    ,  (10) 

then  the  ratio   G',p(v)/  C.-Cv')  depends  only  on  atomic  properties  such  as 
V,  v'  ,  g  ,  etc.,  and  not  on  macroscopic  quantities  like  the  temperature. 
On  account  of  this  fact,  we  may  calculate  this  ratio  of  cross-sections  from  the 
principle  of  detailed  balance  for  conditions  of  statistical  equilibrium,  and 
be  sure  that  the  result  will  still  hold  when  there  is  no  statistical  equilibrium. 
The  condition  of  detailed  balance  is 

^1  ^12  ^'^^•^^^^•^•'^^  "  ^^2  ^21^'^'^*^^^'^*^'''*'''   '  ^-^^ 

where  the  conditions  (?)»  (lO)  must  be  satisfied.  Equation  (lO)  ai-ves  the 
resiilt  v.dv  »  v'.dv'  .  If  we  substitute  a  Maxwellian  function  for  the  velocity 
distribution  function  f(v): 

f(v)  =  Un  (Mj^/2nkT)^/^  ,  v^  ,  exp  (  -  Mj^v^/2kT  |  ,  (12) 
then  equation  (11)  leads  to  the  result 
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(S-^J^^)/  cy^-^{y')   =  (gg/g^).  (v'/v)^  (13) 

*  t 

Substituting  (13)  in  (9),  we  get  *  ' 


g^       ^21^ 


(|tN2)coll     -     ^•^('^)     p-^2     (1^1     ^''       )|      ^  (1^) 

where  ^ 

00 

F(T)  -  /  dv.v.f(v).  0-3^2^"^^   •  ^^^^ 

o 

Neither  the  precise  form  nor  even  the  numerical  value  of  Cr^p(v)  is  known, 

but  in  view  of  the  Boltzmann  factor  in  f(v)   the  precise  fonn  does  not 

much  matter  as  long  as  we  can  make  a  reasonable  estimate  of  CT  close  to  the 

threshold  v  ,  If  we  put 
o 


I  -  !  V     <      V 

0-12  (^)  -  ^  °  (16) 

V    >     V 

o 


then 


F(T)  =  20-  (^)   ^[-e2^/kT]  (1*^).  (l?) 


It  is  clear  that  in  equilibriijm,   substituting  ( 7)  in  (lli)  yields  (^  N^)   ^„-,.,=  0  , 

T 


'8t  2'   coll 
as  must  be  the  case. 


In  writing  down  equation  (lU),  we  asstune  the  form  (12)  for  f(v),  which  implies 
equilibrium  conditions.  Since  we  have  made  the  restriction  (3)  ^md  have 
assumed  the  possibility  of  defining  a  temperature  T  in  (5),  the  use  of 
(12)  in  obtaining  (Hi)  is  a  consistent  assumption. 
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2,2  The  escape  of  energy  by  radiation 


The  rate  of  decrease  of  N_(r,t)  due  to  emission  of  the  quantum 


hv^^  is 


NgCrjt)  I  A^^   +  B^-L  P  (V2-Li£,t)]   -  N^Ujt)   Y(r,t)  ,      (l8) 

where     p     is  the  radiation  energy  density  of  frequency    V     ,  measured  per 
unit  frequency  range.      Instead  of  calculating  this  rate  by  evaluating  the 
density       P     due  to  the  radiation  emitted  by  other  atoms  in  the  system,  we 
proceed  as  follows. 

Let  the  absorption  ciross-section  of  radiation  of  frequency    v     be     k 
so  that  the  radiation  intensities     I  (v),  I  (v),  measured  a  distance     s       apart, 
are  related  by 

Ig(v)     -     I^(v).exp(-Vc^.N^.3)  .  (19) 

ky  is  some  function  of  frequency  v  ,  the  precise  form  depending  on  the 
mechanism  responsible  for  the  line  shape.  Let  P(v)  be  the  intensity  distribu- 
tion of  the  spectral  line,  so  normalized  that 

y  dv.P(v)  a  1,  (20) 

SL 

where  "  SL"  indicates  that  the  integration  extends  over  the  whole  of  the 

spectral  line. 

Now  P  (5,t),  the  probability  that  a  quantum  of  radiation  traverses 

a  distance  ^  =  |  r  -  r'  |  before  being  absorbed,  is 

For  the  fouaulation  of  a  similar  but  less  general  problem,  see  Holstein"^-" 
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r'U,   t)  =  /  d  V  .  P(v).  exp  [-  k^  /  d?.N^(4,  t)  j  .       (21) 

3 
The  number  of  quanta  emitted  in  the  volime  element  d  r'  ,  and  directed 

into  the  solid  angle  dSA  t   per  second,  is 


dS 


^  .  N^Cr'  ,t).  Y  (r',  t).dV   ,  (22) 


and  thus  the  number  of  quanta  that  are  absorbed  in  the  volume  element 
d^^  =  dS.d4     is 

--^  .  N^CrSt).  Y  (rS  t).  d^r'    .  ^A|^  .  d?  .  (23) 

hnK 

The  number  of  quanta  absorbed  per  unit  volume  element  per  second  at   r  due 
to  radiation  coming  from  all  space  is 


d-'C 
■which  is 


■r-  .       /    d  r'      •     expression  (23)   , 


-  (lAn)   .   y    dV    .  N2(r',t)   .     y  (^  »^^*     ^^^H^     *  ^^U) 

G(r,r'jt),  the  probability  that  a  quantum  emitted  at  r'   is  absorbed 
in  tmit  volume  at     r     ,  is  given  by 


G(r,r';t)     =     G(r»,r5t)     -    ^,±£^2^     .  (2^) 

Thus,   from  equations  (18),   (2U),   (25),   the  net  rate  of  change 


of  N„(r,t)   due  to  radiative  processes  is 

(|t  N2)  j.^^  -  -N2(r,t)  Y  (r,t)  +y dV.N2(r' ,t).Y(r',t).G(r,r»;t)  ,     (26) 
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and  combird.ng  this  with  (lii)  we  have  the  total   rate  of  change  of  Np(r,t): 


^\'     X.F(T)[n-N2      (^     *     (gi/gg).^^^       }1      -N2.Y(r,t)  + 


+  y    dV.N2(r',t).Y(r',t).G(r,r'jt) 


which  is  one  of  our  fundamental  equations. 

This  differential-integral  equation  can  be  simplified  somewhat  if 
we  neglect  the  induced  emission  of  radiation  compared  with  the  spontaneous 
emission,   so  that  equation  (l8)   gives 

Y  (r,t)     =     Ap-      "     constant,  (28) 

This  neglect  is  justified  il   there  is  no  incident  radiation  from  an  external 
source,  and  if  p  arising  from  the  atoms  in  the  system  is  small.  If  we  use  (28), 
equation  (2?)  becomes 

N^  -  X,F(T)   [  N  -  N^  I  1  +  ig^/g^).e    '^^        jj  -  N^..^^^  + 


A21  y  dV.N2(r',t).a(r,rJjt) 


(29) 


However,  equation  (29)  still  contains  two  -unknown  functions,  N-(r,t)  and 
T(r,t),  and  thus  to  solve  the  problem  we  need  another  equation,  which  is 
furnished  from  the  definition  (5)  of  the  effective  temperature  T  . 

2.3  Energy  balance  and  the  »  temperature  "  T 

The  kinetic  energy  per  unit  volume  of  the  gas  is  (3/2)kTX,  This 
is  changed  only  by  collisions  and  not  by  radiative  transitions,  because  the 
radiative  processes  by  themselves  cliange  only  N_  and  N^  but  do  not  change  the 
kinetic  energy.  Hence,  the  variations  of  the  temperature  at  any  given  position 
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ajid  time   are  deteimined  by  the  condition 


1^         (3/2)kTX  ]       +         e^^  (|^    N^)     ^^,,     =    0     , 


and  from  the  equation  (ih)  this  reduces  to 

^  e„  At  1 

-  0.   (30) 


1^    [(3/2)kTxj     +     e^i     ^(T)    [    N       -    N^ll     +     (g-L/g2).e     ^^        } 


In  equations  (29)>   (30)  we  have  two  coupled  non-linear  integro- 
differential  equations  which,  together  with  the  appropriate  boundary  conditions, 
will  specify  the  two  functions  N„(r,t)   and  T(r,t), 

3.  Solution  of  a  Simplified  Problem 

The  solution  of  the  two  equations  (29),   (30)  in  their  full  generality 
presents  considerable  difficulties*     In  the  present  work,   as  a  first  approxima- 
tion to  the  solution,  we  neglect  the  effect  of  imprisonment  of  resonance  radia- 
tion by  putting  G(r,r*,t)  «  0,     It  is  clear  that  this  will  give  an  over-estimate 
of  the   rate  of  energy  loss  due  to  radiation  for  the  present  problem  since  in 
this   case  all  quanta  radiated  are  immediately  lost,   there  being  now  no  mechanism 
for  retaining  them  within  the  system. 

To  reduce  the  complications  still  further,  we   tcke  the  total  densities 
N  ,  X  as  absolute  constants,    so  that  there  is  now  no  space  dependence  left  in 
the  problem.     Thus,  our  basic  equations  obtained  from  (29)   and  (30)  are 

I^N^     +     [xF(T)    ^1     +     (g^/g2).e  ^^        J     +     A^-L    1   N^  -  XNF(T)^    (31) 


At 


^    I  (3/2)kTX  }      +   e  21  XF(T)    r    N  -  N^    |l  +   igj^/g^).^  ^^        ] 


0,    (32) 


-  n  - 

where  now  N-  and  T  are  functions  of  tine  only,  and  eII  other  quantities 
(X,  N,  g^,  g^,  e  2i»   ^21^  ^^"^   constants. 

We  now  wish  to  solve  the  problem  of  equations  (31)  and  (32)  for 
t  >  0  with  the  following  initial  conditions 

T  «  T  ,  BT/at  =0,  at  t  =  0  (33) 

In  order  to  solve  the  problem  of  (31)  -  (33),  we  introduce  the 
folloiglng  dimensionless  variables  and  constants: 


n      = 

N^A 

X        = 

e^iAT 

r  = 

c  t 
o 

■)    - 

c^.9^  (x) 

g 

g-^/eg 

11    » 

2N/3X 

\    = 

^2l/^o 

XF(T)       -       c.<^  (x)  (3U) 


where     c       is  a  constant  of  the  order  of   the  reciprocal  of  the  time  between 
o  ^ 

two  excitational  collisions. 

In  this  new  notation,   the  problem  is 

|^+       [    \  W(x).(l  +  g.e^)]     n=9((x)  (35) 

(l/x)     +     \x  ^{x)    [  1  -  n.(l  +  g.e^)  ]  =  0  (36) 

to  be  solved  for  ^    >    0  with  the  initial  conditions 


Substituting  (33)   in  (31)   and  (32),  we  find  that  at  t  =  0,  we  have  N„=N° 
of  equation  (?),   and  also  d  N^/St  =  -  Api^o*     Both  these   results  are 
eminently  reasonable,    and  show  that  the  initial  conditions  (33)  are  satis- 
factory for  describing  how  radiative  2-1  transitions  disturb  the  initial 
statistical  equilibrium. 
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X  -  x^,   8x/3t  -  0,       at  t  -  0,  (37) 

We  shall  take 

(2f  (x)^  e"""  (38) 

as  an  approximation,  bearing  in  mind  that  while  the  precise  form  of 
F(T)  is  not  known,  yet  an  approximation  such  as  (38)  does  yield  the 
basic,  exponential,  T-dependence  of  equation  (17) • 
3.1  Solution  of  the  system  (3$)  -  (37) 

A  standard  way  of  solving  two  coupled  equations  is  to  eliminate 
one  of  the  two  dependent  variables:  here  it  is  clearly  convenient  to 
eliminate  n  rather  that  x  ♦  If  we  carry  out  the  elimination,  we  get 


^  ^  0/3X)(l/x) 
H0  (x) 


1+g.e 


n  (36') 


and 


2       2  /^  ^  2 

dx/dx.,)2     1    I    dx,        -X   LLgx     ,   2  -x   ^ 

-T^^dt^  y^': ^Pdr^*^^^   .M___xuxe   -0. 

dX  I      1+g.e  J  1+ge 

^  (39) 

Thus  the  problem  is  reduced  to  the  solution  of  the  equation  (39) 
with  the  initial  conditions  (37).  The  independent  variable  t  only  ap- 
pears as  d/dT  ,  so  that  it  is  clearly  convenient  to  transform  (39)  so  as 

to  make  dx/dt  (or  some  f'jjiction  of  it)  the  dependent,  and  x  the  inde- 

2 
pendent  variable  }  and  in  particular,  because  the  coefficient  of  (dx/dT) 

in  (39)  differs  from  -(2/x)  by  a  small  quantity  (namely  l/(l+g,e  )),  we 

introduce  the  new  dependent  variable 


y  -  y(x)  -  i^  .  ^   .  (1,0) 

X       * 


When  y  is  substituted  into  (39)  the  latter  will  have  a  y  -term  of  the  form 
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y  /(l+g,e  ),  5ince  this  term  is  snail,  we  may  neglect  it  and  discuss  the 

solution  of 

y  .  dy/dx  +  y  .  J(x)  -  K(x)  =  0 

x^.J(x)  =  X  +  g  +  e""  +  t^^-"""  ^  J  x^.K(x)  =  Xtie"""        (Ul) 

1  +  g.e"" 

with  the  initial  condition 

y  =  0  at  X  =  x  .  (It2) 

This  equation  (hi)  has  been  integrated  numerically  for  the 
following  values  of  the  parameters: 


1    _  _  o    ..     1 

(U3) 


x^  -  ^  ,  g  -  2  ,  n  =  ^  , 


■<-i  ■^rr^  Tn-7 


1  ,   10  ''  ,   10  ^  ,   10  '  . 


Calculations  have  been  made  and  they  show  that  the  neglect  of  the 

2 

y  -temi  does  not  produce  any  significant  changes,  either  in  the 

general  character  of  the  solution  or  in  its  numerical  values,  for 
the  range  of  parameters  used  here» 

I  For  the  theory  and  the  calculation  to  be  valid,  the  condition  (3) 
must  be  satisfied.  As  X  in  (3ii)  is 

A 


the  condition  (3)  in  terms  of  X  becomes 


*-*12 


While  the  value  \  —  1  is  still  permissible  if  O"  »  G"  ,  the  case  X  =  1 
and  Cr  s  io~    cm   in  Table  1  does  not  very  well  satisfy  the  above 
condition,  and  is  hence  included  in  parentheses  in  Table  1  only  for 
comparison. 
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A  detailed  discussion  of  the  solution  of  equation  (hi)   is  given 
in  the  Appendix;  y  rises  very  sharply  just  above  x  (in  fact,  dy/dx  is  oo 
at  X  -  X  )  to  a  maximum  at  x^,  and  then  falls  off  to  zero  as  x  tends  to  in- 
finity; beyond  the  maximum  y(x)  may  be  approximated  by  7)(x),  defined  as 

K|(x)  =  K(x)/J(x).  ihh) 

Once  we  have     y     as  a  function  of     x  ,  we  can  find    T(x)   from 

/x  p 

dx/    X  ,y(x)     ,   and  then  inversion  gives  x  ■  x(T)« 


o 
A  transfonnation  now  gives  the  tempera txire     T     as  a  function  of  time     t   • 

Fig.  1  gives  T(t)   for  T(0)  =  300°K,    C  =     ]D"       cm     and  a  number  of  values 

of  the  Einstein  coefficient  A-,.     Table  1  gives  T  (10  hours)   for  T(o)  =  300 °K 

and  various  values  of    C     and  A„« 

It  is  of  interest  and  importance  to  see  the  extent  to  which  the 

ratio  Nj,/N,   of  the  concentration  of  the  tiro  states  deviates  from  the  Boltz- 

mann  expression  (?)  for  equilibrium  distribution,    ^^  ^s  in  (6)  we  again 

denote  the  equilibrium  value  by  a  superscript  o  ,  we  have  from  (3U) 

N2(t)  n(x) 


N^°(t)  n°(x) 


O/        N  -,//-,.  X 


and  from  (?),  n  (x)  -  l/(l  +  g  e"^).  Hence,  from  (3$) 


N2(t) 


^      1    -    '^|(x)     eVti  -  (g  +  e-^)/(X  +  g  +  e"^)         (16) 
for    X  >  x^  .     Asymptotically,   for  large  x  (i.e.,   t), 

\/^2°    "*    ^/^'^  *  S)       .  (1x9) 

Thus  it  is  seen'th#t  there  will  be  appreciable  deviation  from  the  Boltzmann 
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distribution  if    X     is  not  too  small  compared  with     g     .     Fig.   2  gives  the 
ratio  Np(t)/N2  (t)   for  g  =  2  and  various  values  of  A-,    ,    (T  for  which  \ 
has  the  same  value       1     • 

U.     Application  to   the  Upper  Atmosphere 

We  shall  apply  the  calculations  of  the  preceding  section  for 

the  simplified  model  to   those  regions  of  the  upper  atmosphere  where  there 

3 
are  oxygen  atoms.     Here  we  have  a  triplet  state     P^^  ^    „     with  the  following 

[2I 
energy  difference  and  transition  probabilities  (Pastemack)  "^  -* : 

^P      ^P^  V     »     68  cm"-'-  A     =     1.68  X  :iD~Vsec  . 

o  1 

^P^    ^?^  V     -     158.5  cm"-'-         A     »     8.8     X  10"^/sec  . 

If  we  consider  the  transition     P,  Pp  alone,    and  for  that  region  of  the 

atmosphere  where  the  total  particle  density  X  is    '^    3«10       cm     ,   (corres- 
ponding to   an  altitude  of  about  100  km),    the  appropriate  values  of  the 
parameters  are  g  «  5/3,  (j.  ■a^  l/U,  x     =  0.77,  which  are  comparable  mth  those 
in  (h3).     The  value  of  the  collisional  cross-section  (7   of  equation  (16)   is 
not  known,   and  accordingly  calculations  have  been  carried  out  for  three 
values  of     C   ,   namely  lO"     ,  lO"     ,  lO"       cm     .     The  value  lO"       cm    may 
perhaps  be  a  plausible  one;   corresponding  to  it,   the  value  of  c     of  (3lt)   is 

'~  Ii  X  10     sec.     We  are  pidmarily  interested  in  the   cooling  of  the  gas  over 
night,   30  that  we  want  to  see  how  the  temperature  changes  in  a  time  '^  10  hours. 
From  the  results  in  Table  I,   it  is  seen  that  the  fall  in  temperature  is  about 
30  degrees.     Thus,   at  this  altitude,   the  cooling  by  radiation   overnight  is 
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not  too  great.  This  result  is  a  consequence  of  the  small  values  of  A^^ 
and  the  ratio  N/X  •   For  high  altitudes  where  the  relative  concentration 

of  oxygen  atoms  becomes  greater  (i.e.,  larger  [i     in  equation  (36)),  the 

3 
effect  of  radiative  cooling  due  to  the  transition  in  the   P  of  the  oxygen 

atom  is  somewhat  greater. 

The  above  result  is  obtained  without  considering  the  effect  of 
the  imprisonment  of  the  radiation.  It  is  clear  that  taking  this  into  account 
will  decrease  the  rate  of  cooling.  In  the  particular  problem  of  the  atmosphere 
discussed  above,  this  effect  can  be  expected  to  be  small.  But  in  similar  pro- 
blems where  the  transition  probability  A  is  large,  it  is  necessary  to 
work  with  the  equations  (29)  and  (30), 


* 


In  fact,  for  this  and  lower  altitudes,  a  good  approximation  to  the  result 
of  the  preceding  section  could  have  been  obtained  directly  from  the  ratio 
of  the  energy  radiation  for  unit  volume  of  the  gas  in  10  hours  (calculated 
on  the  asstmption  of  a  Boltzmann  distribution  at  the  initial  temperature) 
and  the  thermal  energy  density  of  the  gas.  For  the  constants  in  (U3), 
this  comes  out  to  be  0,55  per  cent,  leading  to  a  fall  of  16  in  temperature, 
instead  of  the  12  in  Table  I, 
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T/uBLE  I 


Nocturnal  Cooling 
(a)     The  radiative  transition  probability  A^^    for  different  values  of 
the  cross-section    O"  and  the  parameter    X     .      (N,B, !  A-^   ■  ii.lD       sec 

r2i 

is  a  reasonable  value   (Pastemack)  ^ -^) 
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(b)  The  overnight  drop  in  temperature.  The  temperature  listed  as  a 
function  of  C  and  X  is  the  temperature  after  ID  hours  of  radiative 
cooling,  starting  from  an  initial  temperature  T  -  300  K,  See  also  Fig.  1 
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Appendix.  The  numerical  integration  of  equation  (ljl)» 

We  wish  to  integrate  the  equation 
y^     +     J(x).y  -  K(x)      -     0 

2  T/-    ^  ^    ^       ^     -X  ^  [xg«x  (U.) 

X  J(x)    ■     X  +  g  +  e       +  *-a— 


X  .K(x)  -  Xp-.e' 
with  the  initial  condition 


l+g.e'' 


y-0   at   x  =  x^='|.  (U2) 

The  numerical  integration  proper  is  done  in  three  perts  for 
three  different  regions  (i),  (ii),  (iii),  described  below. 

Region  (i):    x  just  above  x  , 

'  o 

At     x  =  x     *y"0,  so  that  sufficiently  close  to  z  ,    (hi) 
o  o 

becomes 

yg     =    K(x„)»  (A.l) 

which  is  integrated  by  quadrature  to  give 

2 

I   =   K(x^).(x  -  x^)      for  y  «  K(x^)/j(x^)o  (A. 2) 

Region  (ii)t   from  region  (i)  up  to  the  maximum  y,  =  y(x ). 

Here  we  simply  carr^r  out  a  straightforward  numerical  integration. 
Region  (iii);  x  >  x_  , 
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Here  dy/dx  is  negative  but  very  small  compared  to  K/y  or  J, 
and  thus  y(x)  is  rather  close  to  a  function 

7J   (x)  =  K(x)/J(x)  ,  (A.3) 

which  is  the  solution  of  an  equation  obtained  from  (Ul)  by  deleting  the 

term  y. dy/dx  ' ' * 

The  functions  y(x)  and  -n  (x)  are  related  by  the  following 

important  relations: 

dy/dx    >   0     if  and  only  if     7  5"^       (•^•^) 
dy/dx    ^0     if  ^^'^   only  if    y  ^  >)  • 

The  actual  numerical  integration  in  range  (iii)  ^rv^^sents  some 
technical  difficulties  because  the  derivative  dy/dx  is  the  small  difference 
of  the  two  relatively  large  quantities  K/y  and  J,  In  trying  to  carry  out 
the  integration  one  observes  empirically  that  errors  build  up  very  rapidly. 


Note  that   d'>;'(x)/dx   <   0,  and  that  its  value  is  small  compared 
with  K/-?^  or  J(x). 

The  purpose  of  the  remainder  of  this  discussion  is  to  estimate  the  maximum 
error  resulting  from  a  replacement  of  y(x)  by  9'i  (x)  for  x  >  x^  o  It 
is  shown  that  for  the  values  of  the  parameters  used  here,  the  maximum 
error  nowhere  exceeds  30  /o ,  and  thus  does  not  affect  the  general  order 
of  magnitude  of  the  results,  which  is  all  that  is  in  question  here.  In 
fact,  we  see  from  Table  II  that  for  \  =  10~  the  error  is  everywhere 
less  than  1  /o. 


-  22  - 

because  the  derivative  is  the  difference  of  two  large  numbers  whose  first 
four  significant  figures  (say)  cancel.  Furthermore,  the  solution  obtained 
numerically  oscillates  with  increasing  amplitude  about  its  mean  value  as 
errors  build  up.  The  standard  rns  thod  of  dealing  with  this  difficulty  is  to 
use  a  large  niimber  of  significant  figures  and  ;.lso  to  carry  out  evaluations 
at  extremely  small  intervals  in  the  independent  variable  x  •  This  method 
is  not  practicable  here,  because  the  calculation  becomes  entirely  too 
lengthy. 

The  procedure  adopted  here  is  to  work  with  extremely  large  intervals, 
h,  in  X  ( '  large"  in  a  sense  to  be  described  below),  end  to  furnish  upper 
and  lower  bounds  for  the  solution  y(x)  for  x  >  x^  ,  We  now  give  the  four 
steps  in  this  procedure, 
l)  The  function  ^   (x)  gives  a  lower  bound  for  y(x)  for  x  >  x^  . 

The  y  -  and  'fJ-   curves  cross  at  x  =  x^,  where  -g^  =  0,  For 
X  >  -x-i  y  >  ^  >  Oj  hence  -^    <    0,   from  equation  {P,h)i   and  no  further 

crossing  of  the  y  -  and  7f-   curves  is  possible  for  x  >  x.,  because  -?*-  <  0, 

2)  The  dependence  of  oscillations  in  the  solution  y(x)  on  the  interval  h« 
Consider  points  (x^,  yjj»nn^  ^^  ^Vl"  V  ^'   ^n+l»  ^  n+1^* 
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We  3sk:  What  is  the  condition  that  must  be  imposed  on  h  in  order  that 

From  (Ul)  we  have 

(d/dx)y    -   K  /y    -    J  ,  (A.6) 

'    n       n  n        n 

so  that 

Then  (A, 5)  will  follow  from  (/.?)  if  the  condition 

'^'•-^    -     ^fl'   <   0  (,.8a) 

is  satisfied}  and  this  is  indeed  the  C£se  if 

h  >   y^  A  .  (A.8b) 

''n  '  n 

Thus  we  see  that  oscillations  in  the  integration  occur  for  h  larger  than 
the  critical  size  of  equation  (A. 8b).  We  see  fron  Table  II  that  this  critical 
size  is  quite  small,  so  that  it  is  not  feasible  to  try  to  avoid  the  oscilla- 
tions by  using  a  sufficiently  small  value  of   h  ,  since  we  have  to  carry 
out  the  integration  over  an  infinite  range  of  x  . 

3)  An  upper  bound  for  y(x)  determined  by  a  point-by-point  method; 
given  ^  ,,  to  find  an  upper  bound  for  y   • 

From  2)  above,  we  infer  that  when  h  is  sufficiently  large,  then 
the  value  of  (y  -  7  ^  ^^-^   increase  as  the  value  of  (  ^^^  -  7^^.2}   increases. 
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Now,  let  the  correct  solution  of  (Ul)  be  y(x).     From  l)  above,  we  know 
that  y      >    "5?      .     Hence  we  may  define  a  solution  z(x)  of  (A,?)   such 
that  z     >^     jZ     -=W_}   and  now  we  know  that 


7^<    y^    <    z^    .  (A.o) 

Further,  we  know  from  (A, 7)  that  z  is  given  by  the  expression 


^n  -  H  7n.l  *  ^-^n    ^  (  7^,1  -  h.j/  -  U.h.K„   l.   (A.7') 


It  must  be  stressed  again  that  (A,9)j   (A»7')  hold  only  for  h  such 
that  equation  (A. 8b)  is  satisfied.     Numerical  values  of  h,  ^ ,   z  for  typical 
x-values  are  listed  in  Table  II,  for  X  =■  1     and    x=  10     , 

1;)     An  upper  bound  for  y(x)  by  a  step-by-step  method  of  integration. 

The  following  method  is  recursive,  rapid,   and  gives  an  adequate 

numerical  upper  bound  for  y(x)  in  the  present  problem.     It  is  a  modified 

stepwise  integration  of  equation  {hi)  with  an  interval  h  larger  than  the 

minimum  set  by  equation  (/,8b)   , 

Consider  two  points  (x^,  y^,  y^)  and  (x^+iJ^n+l'   ^+1^   satisfying 

equation  (A, 5),     Now  replace     y     _   by  y^   ^     ,  where 


^  n+l  •'^n 


(A,10) 


If  h  is  smaller  than  the  critical  value,  this  method  reduces  to  the  usual 


one  given  in  2)  above. 
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li'  y_+(   ^^n+l*  ^^  ^^^   substitute  it,  in  turn,  in  (A. 10),  and  obtain 


^n+2 


J 1/2) 
^n+1 


n+1 


^n.l 
n+1 


(A.IO') 


This  process  may  be  repeated  until  we  reach  some  point  x^  where  jJl^'       "^  V 


m       I  in 


m 


'm+1 


while  y  '   >  ri    ,     At   this  point,  we  continue  as  follows: 


(lA)  .   (1/2)  _  h 
•'m+1     "'m       h 


il/h)    .   (lA)    h 


m 


^m  -  JT72T 


m 


'm+2 


'ra+1 


m+1 


v+1  -  ;iT7iry 

^ra+l 


(A.ll) 


(A.llO 


and  so  on  (eventually  we  introduce  y,  y      ,...)• 

This  procedure  works  quite  well  in  the  present  problem:  it  gives 
and  upper  bound  because  the  absolute  magnitude  of  the  effective  derivative 
decreases  continually  in  this  procedure j  that  is  to  say,  in  (A.IO'),  (A.ll') 
we  always  have 


y(l/2)  ^       ^^  y(*l/2).  y(l/2)_  j^ 

^n+2     ^  n+2'  ^'^^^^^  yn+2     ^n+l    ^ 


and 


m+2     T.m+2' 


but 


(^lA).  J1A)_  h 
V+2     -^ra+l    ? 


^n+1 

"^n+l" 

J 1/2) 

L 

Vl 

\+l 

"^m+1" 

Ji/M 

. 

^m+l 

<V2'   (A.10-) 


<V2-   (^-11") 


Fxirthermore,  in  the  present  problem  this  method  gives  a  rather  smoothly  varying 
function  y  -  y(x). 

From  the  argioments  1)  to  3)  above  it  is  clear  that  the  method  i;)  does 
in  fact  give  an  upper  bound  to  y(x).  In  practice  we  see  from  Table  II  that 
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it  gives  a  better  bound  for  \  •=  1  than  does   the  method  3)>   although  in 
the  case  X  «  lO"'*  the  method  3)   gives  a  better  (i.e.  lower)  bound.     There 
is  of  course  no  inconsistency  in  this. 

y(x)  is  given  as  a  function  of  x  in  Fig,  3  for  X  =  1  and  10" 
Beyond  the  majcLmum,  we  give  a  curve  fitted  smoothly  to   77 (x). 
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Table  II:     solution  of  equation  {hi)   in  region  (iil),  beyond  the  maxjmum. 


X  =>  1 


X 

h 

0.6 

3.6x10"^ 

1.0 

S.QxlD'^ 

2.0 

lD,2yJ.0~^ 

II 

HI 

riU) 

(^)  iV 

.0383 

O.J4 

0.U9 

.0267 

26.8 

0.86 

.OlOli 

9.6 

7.5 

(%) 


X  •»  ID 


-5 


0.6 

1.0 


0.68x00 


-7 


1.28x10' 


-7 


2.0  1.911x10"'^ 


5.30x10""^        3.1 
3.76x10""^  1 

1.50x10""^  1 


.01 
0.1 
0.1 


In  Column  I,  h  is  the  minimum  value  permitted  by  equation  (A. 8b), 

In  Column  II,  z(x)  is  the  upper  bound  for  y(x)  given  by  the  point- 

by-Doint  method  of  the  appendix,  paragraph  3). 

(1/2") 
In  Column  III,  y^  '   (x)  is  the  upper  bound  of  y(x)  given  by  the  stepv/ise 

integration  of  the:  Appendix,  paragraph  U)* 


Note  how  the  errors  in  III  accumulate  with  increasing  x  . 


For  X  =  id"  ,  the  method  of  paragraph  h)   of  the  Appendix  has 
been  refined  by  using  not  ^  dy/dx,  but  instead  (l/p)dy/dx,  where  p  is 
chosen  such  tJiat  y^n   is  as  small  as  possible  but  still  satisfies  the  con- 


dition 


1  _d_ 
p  dx  n 


.(P) 


dx  ^n+1 


which  is  necessary  for  y    to  remain  an  upper  bound  for  y(x)  . 
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